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Abstract. We present hydrodynamic predictions for the charged pion HBT radii 
for a range of initial conditions covering those presumably reached in Pb+Pb col- 
lisions at the LHC. We study central (6 = 0) and semi-central (fe = 7fm) collisions 
and show the expected increase of the HBT radii and their azimuthal oscillations. 
The predicted trends in the oscillation amplitudes reflect a change of the final 
source shape from out-of-plane to in-plane deformation as the initial entropy den- 
sity is increased. 

1. Introduction 

Ideal fluid dynamics has been used successfully to reproduce many aspects of heavy- 
ion collisions at RHIC [1 . However, even though the hydrodynamic model correctly 
describes the hadron spectra at low transverse momenta, it fails to reproduce the 
transverse HBT radii (R s , R a ) at RHIC [21 13 H- It does, however, yield the 
correct normalized oscillation amplitudes for R s and R [H [H [7]. We use it here 
to predict the expected trends for the evolution of the HBT radii at mid-rapidity 
in (A«200) + (t4«200) collisions from RHIC to LHC, including their normalized 
oscillation amplitudes in non-central collisions. These trends may be trustworthy, 
in spite of the model's failure to correctly predict the HBT radii at RHIC. 

2. RHIC -> LHC 

Extending the hydrodynamic code AZHYDRO [5] from RHIC conditions to the LHC 
regime requires a recalculation of the Glauber initial conditions to reflect the achieved 
higher initial energy density and temperature. Some of the expected trends were 
outlined in [9] where the initial temperature was increased to 2 GeV to show the 
possibility of a transition from out-of-plane to in-plane deformation of the freeze- 
out source in non-central Pb+Pb collisions. Such extreme initial temperatures are, 
however, not likely to be reached in Pb+Pb collisions at the LHC. We here present 
predictions for a more realistic range of initial conditions and a better equation of 
state (EoS) than the ideal guark-gluon gas without phase transition used in [9]. 

Hydrodynamics cannot predict its own initial conditions as a function of y^s, 
but it does provide a unique relation between the initial entropy density profile 
and the final hadron multiplicity. The predicted increase of the final charged 
multiplicity from various models ranges from less than twice to more than 4 times the 
multiplicity at RHIC [TQl [11] . We therefore present our results as a function of final 
charged multiplicity, parameterizing it through the initial peak entropy density So at 
thermalization time r in b = collisions. We cover the range from d ^ h (y=b=0) = 680 

(RHIC) to ^(y=&=0)=2040 (LHC). 
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Figure 1. (Color online) Pion HBT radii for central (6=0) Au+Au collisions as 
a function of transverse pair momentum Kt (left) and of initial entropy density 
sq or final charged multiplicity dl J° h (right). 



We use an EoS that transitions from an ideal quark-gluon plasma phase above 
a critical temperature T c to a chemically non-equilibrated hadron resonance gas 
below T c (RappEoS [8]). The final hadron yields are assumed to freeze out directly 
at T c = 164MeV. The hadron momenta are taken to decouple at energy density 
e dcc = 0.075 GeV/fm 3 , T dcc = 100 MeV. Spectra and HBT radii are calculated from the 
emission function obtained via the Cooper- Frye prescription [12j along this decoupling 
surface. For simplicity, we compute the HBT radii from the space-time variances of 
this emission function instead of doing a Gaussian fit to the two-pion correlation 
function, knowing that this overestimates the longitudinal radius Ri by 20-25% [I]. 
A corresponding overall downward correction should thus be applied to all Ri values 
shown below. 

As the initial entropy density so and temperature To increase we reduce the 
thermalization time To, keeping ToTo constant. This yields a reduction from To = 
0.6fm/c for s = H7fm- 3 , ^- = 680 ("RHIC initial conditions") to r = 0.35fm/c 

for s = 602fm" 3 , ^ =2040 ("LHC initial conditions"). 
3. Central collisions 

Figure Q] shows the pion HBT radii for central Au+Au (Pb+Pb) collisions in the 
out-side-long coordinate system [7] as functions of transverse momentum and total 
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Figure 2. (Color online) Pion HBT radii for non-central (6 = 7fm) Au+Au 
collisions as a function of azimuthal angle <J> for different pair momenta Kt, 
for "RHIC" initial conditions (see text). 



charged multiplicity. There are no dramatic changes, neither in magnitude nor in 
if-r-dependence, of the HBT radii as we increase the charged multiplicity by up to a 
factor 3. The largest increase is seen for R s (by ~ 30% at low Kt) while R a (which 
follows R s at Kt = by symmetry) even slightly decreases at large Kt- Ri changes 
hardly at all. For comparison we also performed calculations with a transitionless ideal 
massless gas EoS [9] (not shown graphically). In this case we see the smallest (< 10%) 
increase in R s , about 10-15% increase in R a (at all Kt), and about 25-30% increase 
in Ri . None of these changes will be easy to measure, but one sees that differences 
in the small predicted changes depend on the EoS, in particular on whether or not 
it embodies a phase transition. The main deficiency of hydrodynamic predictions for 
the HBT radii at RHIC (too weak AV-depcndcncc of R s and R a and a ratio R /R s 
much larger than 1) is not likely to be resolved at the LHC unless future LHC data 
completely break with the systematic tendencies observed so far [3] . 



4. Non-central collisions 



One of the strengths of AZHYDRO is simulating anisotropic, non-central collisions. 
For RHIC initial conditions, although the magnitudes of the HBT radii in central 
collisions were not predicted accurately, their normalized oscillation amplitudes at 
small Kt (which measure the source eccentricity at freeze-out [5]) were correctly 
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Figure 3. (Color online) Cuts along the x and y axes through the freeze-out 
surface for (modified) "LHC" initial conditions: — = 2035, so = 1735fm — 3 at 
To = 0.1fm/c (this particular calculation was done for computation of photons 
and dileptons, using a reduced thermalization time to also capture some of the 
pre-equilibrium radiation). 

reproduced [7]. Their extrapolation to LHC initial conditions may therefore have 
predictive power. 

In Figures [2] and [4] we show the azimuthal oscillations of the HBT radii for 
semiperipheral Au+Au collisions at b — 7im, for both "RHIC" and "LHC" initial 
conditions defined above. At low pair momentum Kt, both R s and R show an 
inversion of the sign of the oscillation amplitude between RHIC and LHC. This is 
indicative of a transition from out-of-plane to in-plane deformation of the source at 
freeze-out (S). Figure [3] shows cuts through the decoupling surface (as well as through 
the hypersurfaces indicating the transition from QGP to mixed phase and from mixed 
phase to hadron gas) along the x and y axes (in-plane and out-of-plane, respectively) 
for "LHC" initial conditions: one observes that the source is initially wider in the y- 
direction (out-of-plane), but later becomes larger in the 2-direction (in-plane). With 
an ideal gas EoS this shape transition requires much higher initial entropy densities 
and temperatures [9]. 

Retiere and Lisa [5] introduced the normalized 2 nd order Fourier components 

R lo.s).2 _ R2 (o.s)(°) ~ R to,s)(l) m 

R% R?M + R m ' R lo R?(0) + R?(%) 

of R? Q s ;(<&), shown in Figure [5] as functions of Kt and s . They showed [5] that the 
zero-momentum limit of R 2 s 2 /R 2 s is a direct measure of the spatial source eccentricity 
at freeze-out: e l ° = 2 lim^ T _^ \R^z/R^o) ■ Using this relationship, the lower left panel 
in Figure [S] shows that the freeze-out source eccentricity flips sign between RHIC 
and LHC, and that at LHC the freeze-out source is elongated in the reaction plane 
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Figure 4. (Color online) Same as Figure |21 but for "LHC" initial conditions. 



direction by about as large a factor as it was still out-of-plane elongated in the RHIC 
case. This is different in runs with an ideal massless gas EoS where even with "LHC 
initial conditions" the final freeze-out source is found to be still out-of-plane elongated 
(although just barely so). 

5. Conclusions 

By varying the initial entropy density to control the final charged multiplicity, we used 
the hydrodynamic model to predict trends for the pion HBT radii from Au+Au or 
Pb+Pb collisions as one moves from RHIC to LHC energies. In spite of a documented 
failure of the hydrodynamic model to reproduce the HBT radii measured at RHIC, the 
model has had great success for most other soft-hadron observables, so the predicted 
trends may still be trustworthy. We find very little variation in the HBT radii for 
central collisions, and whatever small differences we see depends sensitively on details 
of the equation of state, in particular whether or not it embodies a quark-hadron 
phase transition. Clear and characteristic changes are predicted for the normalized 
azimuthal oscillation amplitudes of the HBT radii from non-central collisions, indica- 
tive of a qualitative change of the shape of the source at freeze-out which evolves from 
an out-of-plane elongated freeze-out configuration at RHIC to an in-plane elongated 
shape at the LHC. 
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Figure 5. (Color online) Normalized HBT oscillation amplitudes as a function of 
Krp at RHIC and LHC (top) and as function of so for two values of Kx (bottom). 
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